m 


AD-A257  688 


tw.cz- 1  W-o? 

KA-90-39 
January  1991 

l/jofti 


ION  LOSS  BY  COLLISIONAL  UPSCATTERING 


APPROVED  FOR  PL'FL'C  PEf  EASE 
DiSTRlEUTiC:.:  U 


Marlene  Rosenberg 
EMC2 

San  Diego,  CA  92121 


s 


DT1C 

ELECTE I 
N0V23 1992 


CLEARED  f 

Nicholas  A.  Krall 

Krall  Associates  0CT  1  ^  4 

Del  Mar,  CA  9£0M  uiRtciopArt  ros  friedom  of  iww«maijo'- 

and  SECURiTv  REVIEW  (04S9-P«) 

■•-V.RTMtNr  CF  OEFEN'j? 

t  This  work  performed  under  Contract  No.  MDA— 972— 90-C—Q0Q6  for  c  .*  ** 
the  Defenie  Advanced  Research  Projects  Agency,  Defense  Sciences  Office- 


92-30003 

IIISIIRII^ 


REVIEW  OF  THIS  MATERIAL  DOES  NOT  IMPLY 
DEPARTMENT  OF  DEFENSE  INDORSEMENT  OE 
FACTUAL  ACCURACY.  OR  OPINION, 


i®  ENERGY/MATTER  CONVERSION  CORPORATION 
!  9100  A  Center  Street,  Manassas,  VA  22110,  (703)  330-7990 


9  2 


%-s-  wr/ 


1.  INTRODUCTION 


In  at  least  one  idealized  version  of  the  Polywell'VSCIF  (Spherically 
Convergent  Ion  Focus)  fusion  scheme,  the  ions  are  injected  at  a  point  xQ  such 
that  the  magnetic  field  is  sufficiently  reduced  from  its  maximum  value  that  the 
ion  orbits  are  nearly  straight  lines  as  they  converge  to  the  center  of  the 
device.  Typically,  this  location  will  have  a  B-field  in  the  range  of  100  Gauss, 
compared  to  a  Bfflax  of  several  kG.  This  convergence  in  a  strong  potential  (e0 
»  MV^)  produces  a  highly  non-Maxwellian  ion  distribution,  with  a  small  dense 
core  of  nearly  monoenergetic  ions  surrounded  by  a  large  mantle  of  less  dense 
plasma  made  up  of  ions  whose  velocity  is  primarily  radial. 

There  are  two  very  distinct  effects  that  collisions  have  on  ion 
confinement.  Firstly,  as  the  ions  pass  through  the  geometric  center  of  the 
device,  collisional  processes  will  alter  their  velocity  distribution  from  nearly 
monoenergetic  to  Maxwellian.  The  lower  end  of  this  Maxwellian  will  continue  to 
follow  nearly  straight  line  orbits,  but  the  higher  velocity  ions,  making  a  larger 
radial  excursion  into  region  of  higher  magnetic  field,  will  be  deflected  by  the 
B-field  and  will  no  longer  converge  to  the  center.  In  that  sense,  they  are 
"lost"  to  the  dense  core  which  produces  the  bulk  of  the  fusion  reactions  in  the 
Polywell'VSCIF  device. 

There  Is  a  second  collisional  effect.  As  the  ions  traverse  the  mantle 
of  lower  density  plasma  outside  the  core,  collisions  can  transfer  momentum  from 
primarily  radial  flow,  which  produces  good  spherical  convergence,  to  a  more 
Isotropic  distribution  of  velocities.  This  means  an  increase  in  local  azimuthal 
velocity  for  the  ions  outside  the  core,  which  decreases  convergence.  The 
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decrease  In  convergence  Is  a  consequence  of  conservation  of  angular  momentum, 
rv^  »  constant  and  energy  v^  +  v^  -  2e$/M  »  constant,  which  limits  the  minimum 
radius  an  Ion  can  access.  This  effect  will  be  treated  in  a  subsequent  report. 

The  calculation*’^  of  velocity  space  diffusion  due  to  collisions  is 
straightforward;  the  standard  results  from  the  literature  has  been  used  to 
estimate  the  rates  at  which  collisional  processes  in  the  dense  core  and  external 
mantle  produce  diffusion  of  the  ions  In  velocity  .  The  purpose  of  this  note 
is  to  translate  those  calculations  Into  estimates  of  the  loss  rate  of  fuel  ions 
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2.  FORMULATION  OF  THE  EFFECTIVE  LOSS  RATE  DUE  TO  UPSCATTER 


Previous  calculations*  have  shown  that  the  change  in  ion  velocity  in  a 
single  pass  through  the  dense  core  is  small.  Therefore  the  upscattering  should 
be  treated  as  a  diffusion  process  (In  velocity  space)  rather  than  a  single 
scattering  process.  In  analogy  with  the  spatial  diffusion  problem,  we  use  a 
continuity  equation  in  velocity  as  the  basis  for  the  loss  process 
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where  D  -  (Sv)  /r,  with  £v  being  the  change  in  velocity  produced  in  a  time  r  due 
to  scattering  in  the  dense  core,  and  n  the  time  dependent  density  of  ions  in 
velocity  space. 

Density  in  velocity  space  is  the  density  of  interest  to  ion  confinement 
in  the  PolywelltB/SCIF  device,  since  it  determines  the  ion  convergence  and  thus 
the  ion  density  in  the  core.  The  notion  of  ion  loss  in  this  device  is  unique 
among  fusion  schemes;  loss  means  a  decrease  in  ion  core  density  resulting  from 
loss  in  spherical  convergence.  Loss  of  ions  from  the  device  Is  not  a  significant 
energy  drain,  since  source  Ions  have  very  low  energy.  In  fact,  if  an  ion  fails 
to  converge  to  the  core,  It  Is  better  off  lost  from  the  system,  since 
nonconvergent  ions  which  linger  in  the  system  i^itralize  electrons,  reducing  the 
potential  well. 
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3.  CALCULATION  OF  THE  DIFFUSION  COEFFICIENT 


The  diffusion  tensor  in  velocity  space  used  in  the  Fokker-Planck  equation 
Is  given  by 


D 


ij 


It  <Avi4vj> 


(1) 


where  Av  is  the  change  in  the  velocity  of  a  test  particle  in  a  time  r  due  to 

collisions  with  field  particles.  For  an  Isotropic  distribution  of  field 

particles,  the  diffusion  tensor  becomes  diagonal.  It  can  be  expressed  in  terms 

of  components  parallel,  D  ,  and  perpendicular,  Dx,  to  the  direction  of  the 

initial  velocity  (before  the  collision)  of  the  test  particle.  For  reference, 

the  derivation  of  the  components  of  the  diffusion  tensor  when  the  background 

field  particles  can  be  described  by  an  isotropic  Maxwellian  is  given  in  Appendix 

A,  following  Reference  1.  In  the  case  that  the  field  particles  are  ions,  denoted 

2 

by  a  subscript  2,  with  a  temperature  Tz  defined  by  (3/2)T2  -  m1v2/2,  the 
components  of  the  diffusion  tensor  are  given  by 


e2e2n2lnA  ^(Vj/Vg) 
11  ^c^Vj  2(v1/v2)2 

2  2 

e^eJn^nA 
ljp  8«o"lvl 


*  (vj/v2)  - 


2  (v,/v,)4 


(2a) 


(2b) 


4 


Here  *  Is  the  error  function,  defined  by 


*00  -  “  J  exp(-f2)df  ,  (3a) 

ft  o 

and  #j(x)  Is  defined  by 

*i(x)  -  *(x)  -  ^  exp(-x2)  -  *(x)  -  x  ^  .  (3b) 

ft 

Figure  1  shows  a  plot  of  fj(x)/2x2  and  *(x)  -  #j(x)/2x2  (Ref.  2).  In  the  limit 
that  x  - 

*j(x)  -  1 


*(x)  -  1 


For  x  -  1, 


*j(x)/2x2  -  0.214 

*(x)  -  ^(xj^x2  -  0.629  .  ^ 

As  a  first  approximation,  we  assume  that  the  dense  ions  in  the  core  can 
be  described  by  a  Maxwellian  distribution  with  a  temperature  T2  corresponding 
to  the  radial  energy  of  these  energetic  Ions.  We  expect  that  this  estimate  for 


5 


the  diffusion  coefficient  would  be  accurate  to  at  least  about  a  factor  of  2 
compared  with  other  relevant  field  particle  distributions,  such  as  delta 
functions  In  velocity.  For  example,  the  diffusion  coefficient  for  a  test 
particle  interacting  with  an  Isotropic  distribution  of  field  particles,  each 
having  the  same  magnitude  of  velocity,  was  derived  In  Reference  6.  In  a 
coordinate  system  in  which  one  axis  is  parallel  to  the  initial  velocity  vector 
v  of  the  test  particle,  the  diffusion  tensor  D- .*  is  diagonal,  and  is  given  by 

■  J 


where  again  the  subscripts  1,2  refer  to  the  test,  background  particles, 
respectively.  Comparing  the  expressions  for  D„,  Dx  given  in  Eq.  (2)  for  a 
Maxwellian  distribution  of  field  particles  with  the  expressions  for  D„,  Dx  given 
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In  Eq.  (4)  for  a  delta  function  distribution  of  field  particles,  we  obtain  the 
following  estimates.  For  Vj  »  V2,  we  have  approximately  that  D  /D  -*  2/3  and 
DX/0X  -*  1.  For  Vj  -  V£,  we  have  approximately  that  5  /D  -1.6  and  Dj/D^  -  1.1* 

We  note  that  the  use  of  Eq.  (4)  might  be  a  better  approximation  for 
experimental  conditions  in  which  the  core  density  is  in  the  lower  density  range 
of  10*8-l(j^Pcm’8,  while  the  use  of  Eq.  (2)  might  better  apply  to  the  high 
density  case,  n  -  10^-10*8  cm*8. 

3.1  PARALLEL  DIFFUSION  DUE  TO  COLLISIONS  IN  THE  CORE  (UPSCATTER) 

In  this  section  we  consider  the  effects  of  velocity  upscatter  due  to 
scattering  in  the  core.  Upscatter  in  the  radial  direction  due  to  collisions  in 
the  core  would  affect  the  ion  loss  because  an  upscattered  ion  could  make  a  larger 
radial  excursion  into  regions  of  larger  B  field,  where  it  could  be  lost  to  the 
device  in  the  sense  described  above.  We  use  the  parallel  diffusion  component 
D  to  describe  the  diffusion  in  radial  velocity  due  to  collisions  in  the  core. 

The  dense  core  of  the  Polywell”/SCIF  device  occupies  only  a  small 
fraction  of  the  total  volume  of  the  machine.  Since  we  are  only  interested  in 
parallel  diffusion  due  to  collisions  in  the  core,  we  multiply  expression  (2a) 
by  the  fractional  time  that  a  test  Ion  spends  In  the  core  per  pass  through  the 
device.  It  is  assumed  that  the  core  has  radii#  r£  and  uniform  density  ng.  The 
transit  time  through  the  core  is  given  approximately  by  tc  -  2rc/Vj,  where  Vj 
Is  the  velocity  of  the  test  Ion  in  the  core  region.  The  time  for  a  test  Ion  to 
transit  through  the  entire  device,  however,  is  given  approximately  by  tp  ■ 

Vj,  where  xQ  Is  the  radial  location  at  which  the  test  ion  was  born,  and  Vj  Is 
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some  average  velocity  of  the  test  Ion  In  the  device.  Thus  the  effective  parallel 
diffusion  coefficient  due  to  collisions  in  the  core  is  assumed  to  be  given  by 


(5) 


We  rewrite  Eq.  (5)  in  terms  of  the  ion-ion  collision  time  in  the  core, 
defined  by 


Tiic 


e^e^lnA 


(6) 


(Ref.  1).  Using  this  expression,  we  have  that 
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¥*:> 


v2  3  2 


rc  V1 
(rH  v1) 
xo  V1 


1 

T1ic 


(7) 


We  will  use  the  above  expression  to  approximate  the  parallel  diffusion 
coefficient  for  test  Ions  resulting  from  collisions  In  the  core  of  the  device. 

We  need  to  write  the  diffusion  coefficient  In  terms  of  system  parameters. 
To  do  this,  we  use  the  results  of  Ref.  7  for  the  motion  of  an  ion  In  a  slab  model 
of  the  Polywell^/SCIF  device.  From  that  paper,  the  field  particle  Ion  velocity 
Is  given  by 
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vx  '  vxo  •  2e  W  -  #(x0)]/»,  -  -  X-1]2  ,  (8) 


where  x  denotes  the  radial  axis  situated  midway  between  the  magnetic  cusps,  xQ 
Is  the  radial  position  at  which  the  Ions  are  born  with  a  velocity  vXQ,  d(x)  is 
the  electrostatic  potential  which  is  assumed  to  decrease  as  x  decreases,  R  is 
the  radius  of  the  device,  uci  -  eBg/m.,  and  the  magnetic  field  varies  as  B  * 
B0(x/R)m,  where  BQ  is  the  magnetic  field  at  R.  In  order  to  proceed  with  first 
approximation  estimates,  we  assume  that  most  of  the  energy  of  the  core  ions 
arises  from  their  acceleration  by  the  radial  electrostatic  field,  with  the  term 
dp  2ed(x^-  being  the  dominant  term  for  v2  in  the  core.  (It  is  assumed  that 
Ze^x^m.  »  which  is  the  condition  for  ion  reflection  to 

be  avoided  at  the  location  where  the  ions  are  born,  from  Ref.  7.)  Thus  we 
approximate 


2  2e*(x  ) 

v2  *  nfTj 


(9) 


Since  the  test  particle  is  actually  one  of  the  core  particles,  we  assume 
that  Vj  -  v2(l  +  f),  where  f  •  AVj/v2,  and  AVj  is  the  upscatter  In  the  test 
particle  velocity  due  to  collisions  in  the  co#.  Then 


i  [*,(1  +  f)3  — - — 5  1016  (f2)2  i  (10  kV) 
41  (1  +  f)3  s 


r  v,  n,(1012cm"3) 

*  )  — 372 - 

xo  V1  d  '  (10  kV) 


(10) 


9 


for  m.  -  2  nip  for  deuterium.  The  loss  rate  due  to  upscatter  would  then  be  given 
by 


Tloss 


■ 4f2(1  +  f)3  hittW 


<$3/2  (10  kV)  e 
n2(1012  cm*3) 


(11) 


We  estimate  the  value  of  f  which  could  lead  to  ion  "loss,"  in  the  sense 
of  a  broadened  r  ,  due  to  increased  deflection  by  the  higher  B-field  as  the 
upscattered  ion  makes  a  larger  radial  excursion.  Jlo  do  this,  we  estimate  the 
new  radial  location  x'Q  at  which  the  upscattered  ion's  radial  energy  decreases 
to  its  birth  energy.  Neglecting  scattering,  we  assume  the  ion  radial  energy 
profile  is 


E 


r 


Eo  + 


max 


(12) 


where  typically  p  Is  related  to  the  magnetic-field  profile  by  p  -  m.  Due  to 
upscatter  in  the  core,  the  Ion  receives  a  kig^in  energy  AEr, 


xo  p 

<r>  * 


(13) 


Then  E_  +  AE_  -  E.  when 
r  r  o 
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(H) 


The  perpendicular  deflection  of  the  ion  by  the  magnetic  field,  AV,  at  the 
ion  birth  location  is  related  to  the  core  radius  by^ 


rc  AV 

jr  “  r 


(15) 


The  increase  in  core  radius,  Ar  ,  is  then  related  to  the  change  in  AV, 

V* 

AV(x')  -  A V(x0),  by 


Arc  AV(xJ)  -  A V(xQ) 
rc 


(16) 


The  perpendicular  deflection  of  the  Ion  by  the  magnetic  field  was  estimated  in 
Ref.  7,  viz: 


AV(x0) 
-y - 


x _  2m 
cl  ,  o. 


(m  +  1) 


[2e0(xQ)/mi] 


(17) 


At  x' 


I 


AV(x^)  %  j  (1  +  2f)(2m+2/p)x2m+2 

’  ^77?  iP  iw(x007ff,i  •  ( 

so  that  the  change  in  AV  due  to  a  larger  radial  excursion  to  x'Q  is  (f  «  1) 
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Figure  1.  From  Ref.  2,  Fig.  11-6. 
6(x)  •  tj(x)/2x2. 
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